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Abstract. We present the results of spectroscopic analyses of hot horizontal branch (HB) stars in M 13 and 
M 3, which form a famous "second parameter" pair. From the spectra and Stromgren photometry we derived - 
for the first time in M 13 - atmospheric parameters (effective temperature and surface gravity). For stars with 
Stromgren temperatures between 10,000 and 12,000 K we found excellent agreement between the atmospheric 
parameters derived from Stromgren photometry and those derived from Balmer line profile fits. However, for cooler 
stars there is a disagreement in the parameters derived by the two methods, for which we have no satisfactory 
explanation. Stars hotter than 12,000 K show evidence for helium depletion and iron enrichment, both in M 3 
and M 13. Accounting for the iron enrichment substantially improves the agreement with canonical evolutionary 
models, although the derived gravities and masses are still somewhat too low. This remaining discrepancy may 
be an indication that scaled-solar metal-rich model atmospheres do not adequately represent the highly non-solar 
abundance ratios found in blue HB stars affected by diff'usion. We discuss the effects of an enhancement in the 
envelope helium abundance on the atmospheric parameters of the blue HB stars, as might be caused by deep 
mixing on the red giant branch or primordial pollution from an earlier generation of intermediate mass asymptotic 
giant branch stars. 
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1. Introduction 

The two globular clusters M 3 and M 13 form a famous 
"second parameter" pair of clusters that show very dif- 
ferent horizontal branch (HB) morphologies despite quite 
similar metallicities ([Fe/H] = -1.57 for M 3 and -1.54 
for M 13, Harris[r996) : M 3 possesses a horizontal HB pop- 
ulated from the red to the blue end, while M 13 exhibits 
a predominantly blue HB followed by a very long blue 
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tail that extends to temperatures of « 35,000 K (Parise 
ct al. I1998|l . No consensus has yet been reached on the 
reasons for this difference in HB morphology. While it has 
often been argued that age, one of the most widely dis- 
cussed 2"^ parameter candidates (Sarajedini et al. ll997fl . 
cannot be the 2"^^ parameter for this pair (Catelan & de 
Freitas Pacheco ll995l Ferraro et ai-d^SL Paltrinieri et al. 
[TMHl GrundahlUnnnil, the subject remains controversial 
(e.g. Rey et al. I2()()l|l . VandenBerg H2000|l . in particular, 
has emphasized that something other than age must be 
different between M 3 and M 13. 

M 13 (along with oj Cen) harbors the most dramatic 
examples of stars presenting abundance anomalies on the 
upper red giant branch (RGB; Kraft et al. il99^ ITMKl 
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119971 Cavallo & Nagar I2UI)U|I . whereas the abundance 
anomahes in M 3 are considerably less pronounced (Kraft 
et al. [TMl Cavallo & Nagar IMlUjl . In particular, the 
super-oxygen poor stars found near the tip of the RGB 
in M 13 are absent in M 3. These abundance anomalies 
may arise from deep mixing processes which bring up nu- 
clearly processed material from the vicinity of the hydro- 
gen burning shell during a star's RGB phase. Support for 
this possibility comes from the larger Na and Al abun- 
dances and smaller O and Mg abundances found in stars 
near the tip of the RGB in M 13 compared to stars fur- 
ther down the RGB CKraft lT^ Cavallo & Nagar 1^1). 
Alternatively, the recent detections of similar abundance 
anomalies among main-sequence stars in a few clusters 
(Cannon et al. |1998' Gratton et al. '2U(J1^ suggest a pri- 
mordial origin, perhaps due to pollution from an earlier 
generation of intermediate mass asymptotic giant branch 
(AGB) stars. 

Both the deep mixing and primordial scenarios for the 
origin of the abundance anomalies may have consequences 
for the luminosity and morphology of the HB, as outlined 
by Sweigart H1997b|l and D'Antona et al. L2002ji . respec- 
tively. Briefly, in the deep-mixing scenario, helium from 
the H-burning shell is mixed into the envelope on the RGB 
and thereby causes a star to have a bluer and more lu- 
minous position on the HB. In the primordial scenario, 
stars on the main-sequence which have been polluted by 
the products of an earlier AGB generation have a higher 
helium abundance. These polluted stars will then have a 
lower turnoff mass for a given age than the unpolluted 
cluster stars with a normal helium abundance, and thus 
will be more likely to end up on the blue end of the HB. 
As in the deep-mixing scenario, the primordial scenario 
also predicts a higher luminosity of the hydrogen shell on 
the HB, due to the increased helium abundance. 

Thus, both the deep mixing and primordial scenar- 
ios predict that the stars with the strongest abundance 
anomalies should end up on the bluest part of the HB, 
and have both a higher helium abundance and higher lu- 
minosity (and lower gravity) than the HB stars with nor- 
mal abundances. Unfortunately, a straightforward test of 
these predictions is complicated by the processes of dif- 
fusion and radiative levitation in HB stars. The observed 
photospheric helium abundance in hot HB stars, for ex- 
ample, has long been known to be strongly depleted (see 
Moehler i2001i for an overview) presumably due to gravita- 
tional settling. Moreover, the use of luminosity or gravity 
discriminants for testing these scenarios is complicated by 
the supersolar iron abundances in HB stars hotter than 
«11,500 K (Glaspey et al. lH?^ Behr et aipiSI [2000 ) 
due to radiative levitation (Michaud et al. "1983). If the hot 
HB stars with enhanced abundances are analysed with 
model atmospheres at the cluster abundance, the stars 
appear to have anomalously low gravities (Moehler et al. 
I2000II or to be anomalously bright in certain bandpasses, 
such as Stromgren u (Grundahl et al. HM^ . However, 
when Moehler et al. ( 2D0U) used model atmospheres with 
the correct iron abundance for the gravity determinations 



of hot HB stars in NGC 6752, they found that the "low- 
gravity" anomaly mostly disappeared, although a small 
discrepancy remained for stars with 15,000 K < Tes < 
20,000 K. Interestingly, Parise et al. H1998|l also detect a 
possible luminosity offset of M13 HB stars in this temper- 
ature range from their 1620 A photometry obtained with 
the Ultraviolet Imaging Telescope (UIT). 

Although selected hot HB stars in M13 have been 
observed for an abundance analysis (Behr et al. I1999|l 
and studied in Stromgren photometry (Grundahl et al. 
I1999|l . there have been no previous spectroscopic stud- 
ies to obtain temperatures and gravities. Here we present 
spectroscopy of 22 hot HB candidates in M13 to derive 
effective temperatures and surface gravities in order to 
search for any deviations from canonical HB models. We 
also estimate abundances of helium, magnesium, and iron. 
Observations of four hot HB stars in M 3 serve as a control 
sample. 



2. Observations 

We selected our targets in M 13 from the Stromgren pho- 
tometry of Grundahl et al. ( 199H1 see Table and those 
in M 3 from the Johnson photometry of Buonanno et al. 
(1994 see Table We also included three targets in 
M 13 from Behr et al. ifTM??! for comparison. Our tar- 
gets are plotted along the HBs of M 13 and M 3 in Fig. 
1. For our observations we used the Calar Alto 3.5m tele- 
scope with the TWIN spectrograph. While we observed 
in both channels (blue and red) the crowding and conse- 
quent stray light of the primarily red cluster stars made 
the data from the red channel very difficult to reduce and 
analyse. We will therefore limit ourselves to the discus- 
sion of the data from the blue channel. There we used 
CCD#11 (SITe#12a, 800x2000 pixels, (15 /xm)^ pixel 
size, read-out noise 6 e~, conversion factor 1.1 e~ /count) 
and grating T12 (72 Amm^^) to cover a wavelength range 
of 3410 A - 5570 A. Combined with a sht width of 1.5" 
we thus achieved a mean spectral resolution of 3.4 A as 
determined from the FWHM of the wavelength calibra- 
tion lines. The spectra were obtained on June 4-6, 1999. 
For calibration purposes we observed each night ten bias 
frames and ten dome flat-fields with a mean exposure level 
of about 10,000 counts each. Before and after each science 
observation we took HeAr spectra for wavelength calibra- 
tion purposes. We observed dark frames of 3600 and 1800 
sec duration to measure the dark current of the CCD. As 
flux standard stars we used Feige 56, HZ 44, BD-|-28°4211, 
and HZ 21. 

In order to observe as many stars as possible we ori- 
ented the slit to cover up to four hot HB stars at once. 
This of course did not allow us to reduce the light loss due 
to atmospheric dispersion by observing along the parallac- 
tic angle and also required observations in fairly crowded 
regions. 
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Fig. 1. The horizontal branches of M 13 (Grundahl et al. 1998, left panel) and M 3 (Buonanno et al. 1994, right 
panel) with the spectroscopic targets marked by open squares. 



3. Data Reduction 

We first averaged the bias and flat field frames separately 
for each night. The mean bias frames of the second and 
the third night showed the same level of about 923 counts 
(and were therefore averaged), whereas the one from the 
first night had on average 877 counts. The same differ- 
ence was found in the overscan regions of the respective 
bias frames. We therefore determined the mean overscan 
of each science frame before deciding which bias frame to 
use (which was then adjusted to the individual overscan 
level of the science frame). The mean dark current de- 
termined from long dark frames showed no structure and 
turned out to be negligible (1.5±2 counts/hr/pixel). 

We determined the spectral energy distribution of the 
flat field lamp by averaging the mean fiat fields of each 
night along the spatial axis. These one-dimensional "fiat 
field spectra" were then heavily smoothed and used af- 
terwards to normalize the dome flats along the dispersion 
axis. As the normalized flat fields of each night differed 
slightly from each other we always corrected the science 
spectra using the fiat field obtained for the same night. 

For the wavelength calibration we fitted 2"'^-order 
polynomials to the dispersion relations of the HeAr 
spectra (using 46 unblended lines) which resulted in 
mean residuals of <0.15A. We rebinned the frames two- 



dimensionally to constant wavelength steps. Before the 
sky fit the frames were smoothed along the spatial axis to 
erase cosmic ray hits in the background. To determine the 
sky background we had to find regions without any stel- 
lar spectra, which were sometimes not close to the place 
of the object's spectrum. Nevertheless the fiat field cor- 
rection and wavelength calibration turned out to be good 
enough that a linear fit to the spatial distribution of the 
sky light allowed the sky background at the object's posi- 
tion to be reproduced with sufficient accuracy. This means 
in our case that after the fitted sky background was sub- 
tracted from the unsmoothed frame we do not see any 
absorption lines caused by the predominantly red stars of 
the clusters. The sky-subtracted spectra were extracted 
using Home's ( 1986) algorithm as implemented in MIDAS 
(Munich Image Data Analysis System). 

Finally the spectra were corrected for atmospheric 
extinction using the extinction coefficients for La Silla 
fTiig ll977|l as implemented in MIDAS. The data for the 
fiux standard stars were taken from Hamuy et al. H1992I 
Feige 56) and Oke 1)19901 all others) and the response 
curves were fitted by splines or high-order polynomials. 
The fiux-calibration is helpful for the later normalization 
of the spectra as it takes out all large-scale sensitivity vari- 
ations of the instrumental setup. Atmospheric dispersion 
will cause light loss especially at blue end of the spectral 
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Table 1. Coordinates, photometric data, and heliocentric radial velocities for target stars in M 13 (Grundahl ct al. 
1998 [G]; Piotto et al. 2002 [WF3035, WF3485]; this paper [WF3085, WFPC2 data reduced as described by Dolphin 
2000]). H/3 indices with a large (>0™05) uncertainty are marked with a colon. 



Star 


0^2000 


62000 


y 


b — y 


u — b 




[km s^^] 


G43 


16*^41'" 34=7 


+36° 29' 13'.' 7 


14^262 


-0™021 


+0'"688 


2™651 


-287 


G129 


16''41'"41=8 


+36°29'20'.'5 


14"^918 


+0™053 


+1^632 


2™825 


-257 


G137 


16''41'"46=3 


+36° 30' 12'.' 6 


15^019 


+0^040 


+1'^584 


2™837 


-251 


G195 


16''41'"43n 


+36°30'08'.'l 


15^159 


+0™040 


+1'^554 


2™836 


-229 


G201 


16''41'"47=9 


+36°29'47'.'3 


15™237 


+0™025 


+1'^'485 


2™832 


-276 


G208 


16''41"'34n 


+36°29'51'.'4 


15™298 


+0^016 


+1"347 


2™820: 


-244 


G219 


16''41™39n 


+36°29'54'.'0 


15'"222 


+0'"039 


+1'"564 


2™844 


-263 


G235 


16''41'"40n 


+36°29'00'.'l 


15^245 


+0'"033 


+1^528 


2™843 


-245 


G297 


16''41'"35=9 


+36°29'27'.'0 


15^626 


+0'"002 


+1^210 


2™795 


-245 


G314 


16''41'"39=9 


+36°30'00'.'5 


15'"638 


+0'"004 


+1^207 


2™804 


-268 


G322 


16''41'"40=7 


+36°28'52'.'5 


15^865 


-0'"018 


+0'^711 


2^696 


-246 


G342 


16''41'"38=8 


+36°30'42'.'4 


15^840 


-0'"002 


+1'^036 


2™766 


-231 


G382 


16h4i"'4i=i 


+36°29'39'.'l 


16'"058 


-0™021 


+0^649 


2™674 


-255 


G496 


16''41'"4m 


+36°30'05'.'7 


16™538 


-0'"036 


+0^384 


2™673 


-231 


G503 


16''41'''33=8 


+36° 30' 15'.' 5 


16^531 


-0'"035 


+0'"532 


2™713: 


-221 


G635 


16''41™40''5 


+36°31'19'.'7 


16'"957 


-0'"041 


+0™316 


2™696: 


-246 


G744 


16''41'"48=5 


+36°29'59'.'4 


17'"229 


-0'"054 


+0"'266 


2™679: 


-251 


G781 


16''41'"48=8 


+36°29'28'.'2 


17'"348 


-0'"054 


+0"240 


2™661: 


-237 


G827 


16''41'"41=5 


+36°30'37'.'0 


17^530 


-0^072 


+0'^168 


2^671 


-252 


Star 


Q2000 


^2000 


V 


B~V 






^'rad 

[km 


WF3035 


16''41'"48n 


+36°28'20'.'6 


15'5'34 


-0"'01 






-288 


WF3085 


16''41"'37=8 


+36°26'37'.'l 


16'i'15 


-O'^'IO 






-249 


WF3485 


16''41'"36=9 


+36°26'47'.'7 


16'"25 


-0'"06 






-243 



Table 2. Coordinates, photometric data, and heliocentric radial velocities for the target stars in M 3 (Buonanno et 
al. 1994) 



Star 


0^2000 


^2000 


V 


B-V 


frad 












[km S^^] 


B254 


13''42""05"3 


+28°25'09'.' 




-0'^03 


-151 


B352 


13''42™28"8 


+28° 24' 33'.' 


17"01 


-0*^23 


-149 


B445 


13''42™33M 


+28° 23' 57'.' 


16™27 


-0""06 


-146 


B621 


13''41™56"0 


+28° 22' 54'.' 


16';^55 


-O'^^IO 


-177 



range so that the flux distribution of the calibrated spec- 
tra cannot be used to infer temperatures or gravities (e.g. 
from the Balmer jump). 

We derived radial velocities from the positions of the 
Balmer and helium lines. The resulting heliocentric veloc- 
ities are listed in TablesnandEl The error of the velocities 
(as estimated from the scatter of the velocities derived 
from individual lines) is about 30 km s~^. The average 
radial velocity of the M 13 stars of —251 km s~^ agrees 
well with the literature value of —246 km s^^ (Lupton et 
al. I1987P within the error bars. The same is true for M 3 
with —156 km s^^ (our result) vs. —147 km s^^ (Pryor et 
al. |1988). Individual velocities may deviate considerably 
from the mean value because the stars may not have been 
in the center of the slit along the dispersion axis during 
the observation due to the placement of the slit to include 
several target stars simultaneously. 



The Doppler-corrected spectra were then co-added and 
normalized by eye and are plotted in Fig. [3 

4. Atmospheric Parameters 

4.1. Photometric estimates 

To get a first estimate of effective temperatures and sur- 
face gravities we used the Stromgren uvhyfi photometry 
given in Tabled 

Because an empirical calibration of the Stromgren 
indices is not yet available for hot HB stars, we used 
the Stromgren indices (including H/3) computed with the 
ATLAS9 code by Castelli (1998) to derive values for Toff 
and log g. We adopted the model with [Fe/H] = —1.5 with 
a elements enhanced by +0.4 dex, a solar helium abun- 
dance and a microturbulent velocity of 1 km s~^. The 
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normalization of the H/3 indices by Castelli uses the val- 
ues for the Sun and Vega as anchor points for a linear 
interpolation relating the calculated indices, Pcaic to the 
observed indices, Pobs 



/3o, 



bs 



1.287/3, 



calc 



2.5298 



The observed Stromgren uvby photometry was dered- 
dened by Eb-v = 0'^'02. 

For the stars with b~y <0'^01 (wll,000 K), the best-fit 
values of Tcff and log g were determined by minimizing the 

difference between the observed and calculated uvby 



and H/3 photometry. For the stars with b 
computed theoretical (a, r) indices 



y >0 ■ 01, we 



a = 1.36 ■{b~-y) + 0.36 • mO + 0.18 • cO - 0.2448 

r = -0.07 ■{b-y) + 0.35 • cO - (/? - 2.565) 

where for stars with 8500 K < Tcfr< 11,000 K a correlates 
with Tgff and r correlates with logg (Napiwotzki et al. 
I1993|l . The (a,r) method was also used for field HB stars 
with similar Stromgren colours by Kinman et al. (2000). 
The values of Tcff and log g derived from Stromgren pho- 
tometry are given in TableOKCols. 2 and 3) and plotted in 
the upper panel of Fig. 13 These values were then used as 
starting points for the fits of model spectra to the observed 
Balmer lines in order to derive effective temperatures and 
surface gravities (see below). A prior idea of the proba- 
ble range in temperature and gravity is especially helpful 
for stars near the Balmer maximum, where the fit of the 
spectral lines may yield different results depending on the 
starting values, although the solutions on the two sides of 
the Balmer maximum usually show quite different values 
of x^- For reasons of consistency we used the photometri- 
cally determined parameters also as starting values for the 
hotter stars, although the results of the line profile fitting 
depend very little on the starting values for temperatures 
above about 12,000 K. 

For those stars that do not have Stromgren photome- 
try we used B — V from Tables^andl^lto estimate temper- 
atures. As we only have B — V we need additional informa- 
tion to decide on a combination of effective temperature 
and surface gravity. We therefore forced the parameters 
obtained from B — V to agree with the canonical zero-age 
HB model from Sweigart (1997b) for [M/H] = -1.5. These 
parameters (listed in Table O)) should therefore be taken 
with many grains of salt and are useful mainly as starting 
points for the line profile fitting. They are therefore not 
plotted in the top panel of Fig. |3 The B — V colour of 
B 352 in M 3 is very probably wrong, as it suggests an ef- 
fective temperature of well above 25000 K, in which case 
one would expect an absolute brightness of about -|-4™2 
instead of -|-1™9. 



4.2. Line profile fitting 

We first computed model atmospheres using ATLAS9 
(Kurucz il993jl for [Fc/H] -1.5 and [a/Fe] = -f-0.4. 



using a mixing length parameter of 1.25 and no over- 
shoot. This way our models should be identical to those of 
Castelli (1998| used for the analysis of the Stromgren pho- 
tometry and also to those used by Kinman et al. H2000|l 
for the analysis of field horizontal branch stars. Then we 
used Lemke's version^ of the LINFOR program (devel- 
oped originally by Holweger, Steffen, and Steenbock at 
Kiel University) to compute a grid of theoretical spec- 
tra which include the Balmer lines Hq, to H22 and Hei 
lines, but no metal lines. The grid covered the range 
7,000 K < Teff < 15,000 K, 2.5 < logg < 5.5 at a so- 
lar helium abundance log ^^=—1.0. We also computed 
a grid with identical abundances, but no convection, for 
7,000 K < Teff < 35,000 K, 2.5 < log g < 5.5. This grid was 
used for all stars with photometric temperatures above 
12,000 K, where convection should have vanished. 

In Table 121 (Cols. 5 and 6) we list the results obtained 
from fitting the Balmer lines H^ to Hio (excluding H^ 
to avoid the Can H line) with these model spectra. To 
establish the best fit, we used the routines developed by 
Bergeron et al. H1992|l and Saffer et al. H1994fl . as modified 
by Napiwotzki et al. H1999|l . which employ a test. The 
a necessary for the calculation of is estimated from 
the noise in the continuum regions of the spectra. The fit 
program normalizes model spectra and observed spectra 
using the same points for the continuum definition. The 
errors given in Tables|31and^are r.m.s. errors derived from 
the x^ fit (see Moehler et al. il95^ for more details). These 
errors are obtained under the assumption that the only er- 
ror source is statistical noise (derived from the continuum 
of the spectrum). However, errors in the normalization of 
the spectrum, imperfections of flat field/sky background 
correction, variations in the resolution (e.g. due to see- 
ing variations when using a rather large slit width) and 
other effects may produce systematic rather than statis- 
tical errors, which are not well represented by the error 
obtained from the fit routine. As the helium abundance 
is fixed in the model spectra we did not try to fit any 
helium lines, because this would introduce strong biases 
in case of non-solar helium abundances (expected for the 
hot HB stars, see Sect. EJ. The results of this analysis 
are plotted in Fig. 13 (central panel) . As it is not precisely 
clear at which temperature convection truly vanishes in 
HB stars we fitted stars with photometric temperatures 
between 9,000 K and 12,000 K with both the convective 
and the non-convective model spectra. The non-convective 
model spectra yield effective temperatures lower by on av- 
erage 10 K and surface gravities lower by on average 0.015 
dex. Thus the effect of convection at these temperatures 
is certainly negligible for our analysis. 

4.2.1. Stars hotter than 12,000 K 

The star G43 is found to have a high temperature (Teff ~ 
14,000 K) which, along with its bright visual magnitude (y 
= 14™26), suggests that it is a UV-bright star. The star is 



^ For a description see |http://a400. sternwarte.uni-erlangen.de/~ai26/linfit/ 
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Table 3. Atmospheric parameters for our programme stars in M 3 (B) and M 13 (G, WF) as derived from photometry 
(Cols. 2 and 3) and fits to the Balmer line profiles using metal-poor, a-enhanced model atmospheres with solar helium 
abundances. Cols. 5 and 6 list the results obtained with model spectra which contain only hydrogen and helium lines, 
while Cols. 7 and 8 give the results obtained with model spectra that also include metal lines for [Fe/H] = — 1.5 and 
[a/Fe] = +0.4. For those stars where we could obtain a meaningful fit we also obtained iron abundances from spectrum 
synthesis (Cols. 4 and 9, see text for details). Given are abundances using microturbulent velocities of 3 km s~^ and 
(in brackets) for km s~^. The errors are statistical errors only. 



Star Photometry log epc Line profiles log epc 





Teff 

[K] 


lo 

[cm s 


gg 




Teff 

[K] 




lo 

[cm s 


gg 


Teff 

[K] 




lo 
[cm s 


gg 










G43 


13050 


2 


60 




13900±710 


3 


14±0 


14 


14400±450 


3 


24±0 


09 








- 


G129 


8630 


3 


01 


5.7 (5.6) 


8000± 70 


2 


88±0 


09 


7900± 50 


2 


93±0 


07 


5 


1 


(5 


1) 


G137 


8910 


3 


14 


6.4 (6.8) 


8400± 50 


3 


04±0 


02 


8400± 50 


3 


08±0 


02 
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(6 


4) 


G195 


8980 


3 


19 


5.8 (5.8) 


8200± 40 


2 


95±0 


03 


8300± 40 


2 


99±0 


02 


5 


2 


(5 


2) 


G201 


9360 


3 


26 


6.4 (6.6) 


8900±200 


3 


12±0 


12 


8300± 50 


2 


82±0 


03 


5 


5 


(5 


6) 


G208 


9870 


3 


43 


6.6 (6.6) 


10900±200 


3 


73±0 


07 


10900±200 


3 


76±0 


07 


6 


9 


(7 


0) 


G219 


8970 


3 


22 


5.8 (6.0) 


10300±160 


3 


89±0 


07 


8100± 60 


2 


85±0 


05 


4 


9 


(5 


2) 


G235 


9150 


3 


27 


5.9 (5.8) 


11300±130 


4 


12±0 


05 


11300±130 


4 


15±0 


05 


6 


8 


(6 


7) 


G297 


10590 


3 


59 




10600±130 


3 


60±0 


05 


10600±120 


3 


63±0 


05 


6 


5 


(6 


6) 


GSM 


10620 


3 


66 




10400±190 


3 


64±0 


07 


10500±180 


3 


68±0 


07 


6 





(6 


0) 


G322 


13750 


3 


58 




15000±370 


3 


62±0 


07 


14900±360 


3 


61±0 


07 


7 


4 


(7 


5) 


G342 


11520 


3 


69 




11300±140 


3 


71±0 


05 


11300±140 


3 


75±0 


05 


6 


2 


(6 


1) 


G382 


14160 


3 


38 




14300±240 


3 


72±0 


05 


14300±240 


3 


73±0 


05 


8 





(8 


0) 


G496 


17560 


4 


03 




18600±350 


4 


20±0 


05 


















G503 


15800 


4 


31 




14900±580 


3 


86±0 


10 


14900±570 


3 


86±0 


10 


7 


4 


(7 


6) 


G635 


18800 


4 


82 




19500±450 


4 


38±0 


07 


















G744 


19880 


4 


81 




21000±620 


4 


52±0 


10 


















G781 


20390 


5 


00 




19700±820 


4 


61±0 


14 


















G827 


21720 


4 


65 




21200±430 


4 


73±0 


07 


















WF3035 


9550 


3 


51 




10300±240 


3 


74±0 


10 


10300±240 


3 


77±0 


10 










WF3085 


10260 


3 


66 




12900±300 


3 


36±0 


07 


12900±300 


3 


41±0 


07 


7 


5 


(8 


2) 


WF3485 


9790 


3 


57 




13300±330 


3 


76±0 


07 


13400±330 


3 


78±0 


07 


7 


9 


(8 


4) 


B254 


9630 


3 


53 




10600±110 


3 


62±0 


05 


10600±110 


3 


67±0 


03 


6 


3 


(6 


2) 


B352 


28600 


5 


63 




13900±300 


3 


66±0 


05 


13900±300 


3 


69±0 


05 


8 


1 


(8 


6) 


B445 


9850 


3 


58 




9700±170 


3 


28±0 


09 


9800±170 


3 


35±0 


09 


6 


2 


(6 


2) 


B621 


10170 


3 


64 




12900±260 


3 


71±0 


07 


13000±250 


3 


75±0 


07 


7 


6 


(7 


9) 



clearly seen as such on the 1600 A image of M13 obtained 
with the Ultraviolet Imaging Telescope (UIT, Parise et al. 
I1998|l , but the crowding prevents a good determination of 
its ultraviolet magnitude. Using the y magnitude and the 
derived effective temperature, and assuming an apparent 
distance modulus of {m—M)v — 14™51 toward M 13 and a 
reddening of Eb_v = 0™02, we estimate a luminosity of log 
L/ L0= 2.47. The evolutionary status of G43 is uncertain, 
but given the large extreme HB (EHB) population in M13, 
it may be either a He-shell burning post-EHB star, or a 
post-early AGB star (cf. Moehler et al. 1998). In any case, 
since G43 is almost certainly not an HB star, it is excluded 
from further discussion in this paper. 

The results from the Balmcr line profile analyses place 
the stars hotter than 12,000 K further away from the 
ZAHB than the results from the Stromgren colours (see 
Fig. 121 central panel). These stars, however, show also 
weaker He i lines than predicted by model atmospheres 
with solar helium abundance, while in the cooler stars 
the He I lines - where detectable - are well described by 
line profiles for solar helium abundance. The offset from 



the canonical tracks seen in Fig. |3 for the stars hotter 
than about 12,000 K in both M 3 and M 13 is aheady 
well known from other globular clusters (see, e.g. Moehler 
12001(1 . To better illustrate this general behavior, we also 
give the results for the blue HB stars in NGC 6752 ob- 
tained with similar methods by Moehler et al. ( 200^. Note 
that the gravities of the blue HB stars in M 3, M 13 and 
NGC 6752 all follow the same trend with Toff. Possible 
reasons for this behaviour are discussed in Sects. El and El 

4.2.2. Stars cooler than 12,000 K 

As can be seen from the top and central panel in Fig. (21 
and from Table (21 there is an excellent agreement be- 
tween the values of Toft a-nd log g derived from Stromgren 
photometry, and those derived from a detailed Balmer 
line fit for the three stars (G297, G314, G342) in M13 
with Stromgren temperatures between 10,000 K and 
12,000 K. This agreement provides evidence in support 
of the Stromgren calibration of Castelli (IU98') for deriv- 
ing the atmospheric parameters of low-metallicity hot HB 
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Fig. 3. Atmospheric parameters of the programme stars 
in M 3 and M 13 obtained with a-enhanced metal-poor 
model atmospheres ([Fe/H] = -1.5, [a/Fe] +0.4). The 
sohd Unes mark the zero-age (ZAHB) and terminal-age 
(TAHB) loci of canonical HB tracks for [M/H] = -1.54. 
These loci define the region within which the HB models 
spend 99% of their HB lifetime. 

a) Temperatures and gravities from Stromgren uvbyP 
photometry (for comparison we also show results for 
cool field BHB stars obtained with the same method by 
Kinman et al.|2000|. 

b) Temperatures and gravities for all stars are derived 
from fits to the Balmer line profiles, using as starting val- 
ues the photometrically determined parameters. The error 
bars note the formal 1 a errors from the fine profile fits. 
The model spectra used here do not include metal lines. 
For comparison we also show the results for blue tail stars 
in NGC 6752 from Moehler et al. ijMM obtained with 
metal-poor model spectra without metal lines). 

c) Temperatures and gravities of the programme stars 
cooler than 15,000 K from Balmer line profile fits using 
model spectra that include metal lines for [Fe/H] = — 1.5 
and [a/Fe] = +0A. 



stars in this temperature range. The atmospheric param- 
eters of these stars as well as B254 also agree well with 
canonical zero-age HB (ZAHB). 

As can be seen in the top panel of Fig. |2| the stars 
with Stromgren temperatures cooler than R:!lO,000 K show 



a small offset towards low gravities from the canonical 
HB. Fig.^shows that at least the two coolest stars might 
actuaUy be somewhat evolved as they are brighter than 
the main HB in that region of colour-magnitude diagram. 
There is also good agreement with gravities derived for 
field blue HB stars in this temperature range by Kinman 
et al. ifiroO;) . 

Comparing the top and central panel of Fig. |3 how- 
ever, shows that there are significant discrepancies be- 
tween the results obtained from photometric and from 
spectroscopic data for these stars (cf. Table . The re- 
sults from the Balmer line profile fitting move the cool 
BHB stars to hotter and cooler temperatures, thereby pro- 
ducing a gap between about 9000 K to 10,000 K. 

While these stars are not the central target of our in- 
vestigation we followed the suggestions of the referee to 
look for causes that might explain this behaviour. Readers 
interested in the hot stars in M 3 and M 13 may skip the 
remainder of this section and move directly to Sect.El 

The Balmer line profile fitting routine could not find a 
a model spectrum that provides a good fit to all Balmer 
lines in two cases (G235 and WF3035, both below the 
ZAHB in Fig.EJ: Models at about 10,000 K to 11,000 K 
predict too weak H^ to H^ lines, whereas models at about 
8000 K predict too weak Hio to Hg fines. The spectra of 
ah other cool stars are very well reproduced by the model 
fits. We checked the fohowing possible causes for problems 
in our Balmer line fitting: 

Synthetic Spectra: We compared "our" temperature- 
pressure distribution and resulting theoretical Balmer 
profiles for selected models with those computed 
by Castelh (available on http://kurucz.harvard.edu/) 
and those computed using the program SYNSPEC 
(Hubeny & Lanz I2()()3|l and found no discernible dif- 
ferences. 

Metal lines: The Balmer line fitting procedure we used 
was originally developed for hot white dwarfs, which 
show only hydrogen and/or hehum lines. Blue HB stars 
below 10,000 K, however, may show strong metal lines 
even in case of metal-poor stars, which may affect the 
shape of the Balmer lines especially when observed 
at low resolution. To rule out the influence of metal 
lines on the line profile fitting for the cooler stars we 
produced a grid of model spectra with [Fc/H] = — 1.5 
and [a/Fe] — +0.4 (with and without convection for 
the same temperature ranges as above), which include 
metal lines (up to effective temperatures of 15,000 K) 
and which we used for stars with photometric tem- 
perature estimates below 15,000 K. The results of fits 
with these model atmospheres are plotted in the bot- 
tom panel of Fig. 13 and listed in Table 01 (Cols. 7 
and 8). Two stars show large changes (G201 moves 
from 8900 K to 8300 K, G219 moves from 10300 K 
to 8100 K) and the of all stars is significantly re- 
duced. However, the discrepancies with the evolution- 
ary tracks persist as well as the bad fits for G 235 and 
WF 3035. 
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Table 5. Magnesium abundances for our programme stars 
in M 3 (B) and M 13 (G, WF) for microturbulent veloc- 
ities of 3 km s~^ and (in brackets) km s~^. We also 
give the atmospheric parameters used for the abundance 
determinations. 



Star 


T^ff 




loe ^ 








[cm s~^] 


3 km/s (0 km/s) 


G297 


10600 


3.63 


— 1.00 


6.2 (6.5) 


VjrO 




o.uo 


— 1 no 




G322 






—2.32 




G342 


11300 


3.75 


-1.00 


5.9 (6.2) 


G382 


14100 


3.94 


-2.07 


6.1 (6.3) 


G496 


17700 


4.33 


-1.90 


6.2 (6.3) 


G503 


15200 


4.20 


-2.60 


6.7 (6.9) 


G635 


19200 


4.58 


-2.15 


6.9 (7.0) 


WF3085 


13100 


3.65 


-2.34 


6.3 (6.6) 


WF3485 


13300 


4.05 


-2.31 


6.0 (6.1) 


B254 


10600 


3.67 


-1.00 


6.4 (6.7) 


B352 


13800 


3.92 


-2.21 


6.1 (6.2) 


B445 


9800 


3.35 


-1.00 


6.0 (6.3) 


B621 


13000 


3.98 


-2.20 


6.2 (6.4) 



Signal-to-Noise and/or resolution: Our data have S/N 
between 40 and 60 (determined from the continuum 
between 4500 A and 4600 A) and a resolution of 3.4A. 
We convolved a theoretical spectrum from our grid 
with metal lines to this resolution and added noise to 
achieve a S/N of 40 and 60, respectively. Fitting this 
spectrum reproduced the original parameters (9000 K, 
logg= 3.0) to within ±100K in Toff and ±0.02 dex in 
logg. 

Assuming that the real resolution might be better than 
determined from the wavelength calibration frames 
(which illuminate the whole slit of I'.'S) we assumed a 
resolution of 3.0A for another test, both with the the- 
oretical model spectrum as with our data. Again the 
effect is small, with changes of less than 50 K in Tcff 
and 0.02 dex in logg. We also verified that the shape 
of the wavelength calibration lines is well described by 
a Gaussian. 

Rotation: Assuming a projected rotational velocity of 20 
km s~^ does not affect the results, which is not unex- 
pected considering the limited resolution of our spec- 
tra. 

So we currently have no explanation for the disagree- 
ment between the atmospheric parameters for the cool 
stars derived from Stromgren photometry, and those de- 
rived from Balmer-line fitting. 

5. Radiative levitation of heavy elements 

As described in Moehler et al. H2000|l and Behr et al. 
H1999lOr)UU|l . hot HB stars show helium depletion and iron 
enrichment for effective temperatures above 11,000 K to 
12,000 K. To verify this behaviour from our low-resolution 
spectra, we performed a spectrum synthesis to reproduce 
the iron lines in the wavelength region 4490 - 4590 A (as 



3 

N 

6 

u 
o 
G 



4500 4520 4540 4560 4580 

x[A] 

Fig. 4. Normalized spectra of the programme stars in 
M 13 and M 3 and theoretical spectra for the atmospheric 
parameters (from the Balmer line profiles) and iron abun- 
dances listed in TableOl For stars hotter than 12,000 K we 
show the fits obtained for atmospheric parameters from 
metal-rich model atmosphere fits. The solid and dotted 
lines mark the strongest and moderately strong iron lines 
in the spectra. 

described in Moehler et al. I2()()[)|l . First estimates are ob- 
tained using the atmospheric parameters given in Cols. 7 
and 8 of Table and are listed in column 9 of the same 
table. The best fitting theoretical spectra are shown in 
Fig. ^ together with the observed spectra. While it is 
obvious from Fig. 0] that the cool stars show significant 
iron lines our fits (using the atmospheric parameters from 
the Balmer line profile fits) indicate very low iron abun- 
dances. If wc use the atmospheric parameters derived from 
Stromgren photometry instead the iron abundances are 
much closer to [Fe/H] « —1.5 (see column 4 of Table OJ. 
We take this as further evidence that the results from the 
Balmer line profile fits cannot be trusted for stars cooler 
than about 9500 K and we will omit these stars (including 
WF 3035 due to its fit problems) from all further discus- 
sion, as they are not the central goal of this paper. 

In order to be at least roughly consistent with the ac- 
tual stellar abundances, we redetermined Toff, logg and 
log for the stars showing evidence for enriched iron 
abundances using model atmospheres with super-solar 
metallicity ([M/H] — -1-0.5), varying helium abundance 
and no convection. This was done by fitting the Balmer 
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Table 4. Atmospheric parameters for our programme stars in M 3 (B) and M 13 (G, WF) as derived from fits to 
Balmer and helium lines using enriched scaled-solar model spectra ([M/H] = +0.5). Cols. 2-4 give the results obtained 
with model spectra containing only hydrogen and helium lines, Cols. 6-8 give the results obtained with model spectra 
including metal lines. For those stars where we could obtain a meaningful fit we also obtained iron abundances from 
spectrum synthesis (Col. 5, see text for details). Given are abundances using microturbulent velocities of 3 km s^^ 
and (in brackets) for km s^^. The errors are statistical errors only. 



Star 


Line pr 


ofiles w/o metal lines 


lo 


g epo 


Line profiles w metal lines 






Tcff 


logg 








Teff 


logg 


log^ 






[K] 


r — 2i 

[cm S ] 






[K] 


r ~2i 

[cm S ] 


G322 


14600± 


400 


3.82±0.09 


-2.32±0.23 


7.4 


(7.5) 


14500±400 


3.84±0.09 


-2.41±0.28 


G382 


14100± 


200 


3.94±0.05 


-2.07±0.12 


7.8 


(8.2) 


14100±300 


4.00±0.07 


-2.13±0.14 


G496 


17700dz 


400 


4.33±0.07 


-1.90±0.09 






17300±400 


4.35±0.07 


-1.92±0.09 


G503 


15200± 


600 


4.20±0.12 


-2.60±0.42 


7.4 


(7.6) 


15100±500 


4.22±0.12 


-2.72±0.48 


G635 


19200± 


600 


4.58±0.09 


-2.15±0.10 






18700±600 


4.59±0.09 


-2.16±0.12 


G744 


20200±1100 


4.61±0.12 


-2.04±0.12 






18400±700 


4.57±0.12 


-1.94±0.12 


G781 


18800±1000 


4.71±0.16 


-1.46±0.16 






17800±800 


4.70±0.16 


-1.37±0.17 


G827 


20800± 


700 


4.89±0.09 


-2.37±0.10 






19800±600 


4.89±0.09 


-2.28±0.12 


WF3085 


13100± 


300 


3.65±0.09 


-2.34±0.35 


7.4 


(7.9) 


12800±300 


3.62±0.09 


-2.36±0.38 


WF3485 


13300± 


300 


4.05±0.09 


-2.31±0.38 


7.7 


(8.2) 


13300±300 


4.08±0.09 


-2.36±0.45 


B352 


13800± 


300 


3.92±0.07 


-2.21±0.23 


8.1 


(8.5) 


13800±300 


3.99±0.07 


-2.29±0.29 


B621 


13000± 


300 


3.98±0.07 


-2.20±0.31 


7.4 (7.7) 


12800±300 


3.98±0.09 


-2.34±0.36 



lines to Hio (excluding He to avoid the Can H line) 
and the Hei hues 4026 A, 4388 A, 4471 A, and 4921 A 
with these metal-rich model atmospheres. The results are 
given in Cols. 2-4 of Tabled) Using these new atmospheric 
parameters, we then repeated the iron abundance analysis 
and used those fits for Fig. [S| 

We also determined magnesium abundances where 
possible from the equivalent width of the Mgii line at 
4482A (see Table EJ. 

The final iron abundances derived from our spectra are 
hsted in Col. 9 of Table El for stars cooler than 12,000 K 
and in Col. 5 of Table^lfor stars hotter than 12,000 K and 
are plotted in the central panel of Fig. [S| For consistency 
we always plot the abundances obtained with a microtur- 
bulent velocity of 3 km s^^, although km s^^ would 
be more appropriate for the stars affected by diffusion. 
A drastic change in the iron abundance between 11,500 
K and 13,000 K is obvious, in good agreement with the 
findings of Glaspey et al. fl9891 for two hot HB stars in 
NGC 6752 and Behr et al. (1999, 2000) for hot HB stars 
in M 13 and M 15. The results in Fig. El show that aU 
stars hotter than 12,000 K for which we can estimate the 
iron abundance show evidence for radiative levitation. The 
iron abundance for the hotter stars is a factor of 50-100 
greater than that of the cluster in general and consistent 
with that required to explain the Stromgren u-jump dis- 
cussed by Grundahl et al. H1999I logepe = 8.1). We also 
find a one-to-one correspondence between the position of 
a star relative to the u-jump (blueward/redward) and its 
iron abundance (enhanced/cluster abundance) whenever 
we were able to determine an iron abundance (cf. Fig.Ej). 

The upper panel of Fig. El shows that the onset of ra- 
diative levitation is accompanied by a drop in the helium 
abundance by 1 dex or more. As the He i lines are very 
weak at about 10,000 K we did not try to determine he- 



lium abundances for the cool stars. We note, however, that 
the observed He i lines in these cool stars agree well with 
those predicted for solar helium abundance and thus with 
the results of Behr et al. I|1999[ 12000(1 from high resolu- 
tion spectroscopy. Clearly the stable stellar atmosphere 
required for radiative levitation also permits the gravita- 
tional settling of helium. 

The mean magnesium abundance for the stars cooler 
than 12,000 K is [Mg/H] = -1.4 and for the hotter stars 
it is —1.2. So the magnesium abundance does not change 
significantly during the onset of diffusion. The same result 
was also found by Moehler et al. l|^nnn)l for NGC 6752 and 
Behr et al. ifMHH IMMljl for M 13 and M 15. 

Rey et al. ()2001ll have suggested that radiative levita- 
tion and helium diffusion may be responsible for the for- 
mation of a blue tail in M 13. However, we find that the 
HB stars in M 3 show the same abundance pattern with 
respect to effective temperature as the stars in M 13 and 
NGC 6752. The fact that M 3 contains stars affected by 
radiative levitation but does not have a blue tail provides 
strong evidence against the Rey et al. H2001|l suggestion. 

The stellar parameters derived for the stars hotter 
than 12,000 K from the metal-rich model atmospheres 
are plotted in the upper panel of Fig. and compared 
to the canonical HB locus for a helium abundance Y of 
0.23 and a scaled-solar metallicity [M/H] of —1.54. As al- 
ready noted by Moehler et al. H2000|l . using scaled-solar 
metal-rich model atmospheres moves the stars hotter than 
12,000 K closer to the canonical ZAHB, substantially re- 
ducing the low gravity offset seen in Fig. El However, it is 
obvious that discrepancies still exist between atmospheric 
parameters derived from observations and the predictions 
of canonical HB theory. 

The model spectra calculated from the metal-rich 
model atmospheres and used for the results in the up- 
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Fig. 5. Abundances of a) helium, b) iron, and c) magne- 
sium for the programme stars in M 3 (open squares) and 
M 13 (filled squares). Also given are the results of Bchr 
(priv. comm.) for stars in M 3 (open triangles) and M 13 
(filled triangles). 



per panel of Fig.[7|did not include metal lines. In order to 
determine if the strong metal lines that show up at the on- 
set of diffusion might influence our results, we calculated a 
new grid of metal-rich model spectra, which include metal 
lines, for temperatures up to 23,000 K, and then used these 
spectra to fit all stars with effective temperatures above 
12,000 K. The results are plotted in Fig. |7| (lower panel) 
and fisted in Table H (Cols. 6-8). While the differences 
for stars cooler than 16000 K are small, the hotter stars 
become considerably cooler when fitted with theoretical 
spectra that include metal lines. This shift in tempera- 
ture yields consistency of the hottest stars with canoni- 
cal HB evolution tracks. One should remember, however, 
that these spectra were computed using scaled-solar abun- 
dances. It is very unlikely that diffusion will produce an 
overall enrichment of elements - highly non-solar abun- 
dance ratios are much more probable. So these results, 
while promising, should be taken with a grain of salt. 

Also, as can be seen from Fig.|51 the masses of the stars 
obtained from their values of Toff and log g in Tables 
(Cols. 7 and 8) and2](Cols. 6 and 7) tend to be too low. 
Comparing the derived masses to those expected for a star 
on the ZAHB at the same temperature yields a mean mass 



15 - M 13 



16 - 



IB - 



19 r . I • • - i : 

-0,5 0,5 1 1,5 

u-y 

Fig. 6. The w-jump in M 13 with the spectroscopic targets 
marked: Open triangles mark stars with low iron abun- 
dance ([Fe/H] w —1.5), open squares mark stars with high 
iron abundance ([Fe/H] w -1-0.5). Stars for which the iron 
abundance could not be determined due to the limited 
S/N (hot stars) or uncertainties in the atmospheric pa- 
rameters (cool stars) are marked by open circles. 



ratio of 0.9±0.2 for the 5 stars below 12,000 K (analysed 
with metal-poor model spectra including metal lines) and 
0.8±0.2 for the 12 hotter stars when analysed with metal- 
rich model spectra including metal lines. We used distance 
moduh of (m - M)v = U^51 (M 13) and 15'°02 (M 3) 
to derive these masses. 

We conclude that both the gravities and masses of the 
HB stars between 12,000 K and 16,000 K remain about 
0.15-0.2 dex too low compared to canonical predictions 
even when the stars are analyzed with metal-rich model 
spectra which include metal lines. This remaining offset 
may reflect a systematic error caused by our use of metal- 
rich model atmospheres with scaled-solar abundances. The 
results of Behr et al. H1999[Ol)00|l have demonstrated that 
radiative levitation can lead to severely nonsolar abun- 
dance ratios in blue HB stars. Quite possibly, the use of 
scaled-solar metal-rich model atmospheres may not suf- 
ficiently approximate the actual atmospheric abundances 
in these stars. In addition, at such a high ([Fe/H] = +0.5) 
metal abundance, the model atmospheres are not well 
tested, and are more sensitive to inadequacies in the opac- 
ity distribution function than are models at lower abun- 
dances. Also a stratification of the stellar atmosphere in- 
fiuences the line profiles and our model atmospheres as- 
sume homogeneous atmospheres. Another possibility, sug- 
gested by Vink & Cassisi II2002|I . is that an enhanced stel- 
lar wind in the radiatively levitated HB stars may alter the 
wings of the Balmer line profiles, leading to an underes- 
timate of the surface gravity. While it remains to be seen 
if radiative levitation can be effective at the high mass 
loss rates obtained by Vink & Cassisi H2002|l . their results 
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Fig. 7. Temperatures and gravities of the programme 
stars that show evidence for iron enrichment. The sohd 
hnes mark the zero-age (ZAHB) and terminal-age (TAHB) 
loci of canonical HB tracks for [M/H] = —1.54. These loci 
define the region within which the HB models spend 99% 
of their HB lifetime. 

a) Here we adopted a super-solar metallicity ([M/H] — 
-1-0.5) for the model atmospheres (see Sect.|5lfor details), 
but did not include metal lines in the theoretical spectra. 
For comparison we also show the results for blue tail stars 
in NGC 6752 from Moehler et al. ijMjUl obtained with 
metal-rich model spectra without metal lines). 

b) These are the results obtained with metal-rich model 
spectra that include metal lines. 

do raise potential concerns about the use of hydrostatic 
model atmospheres. 

6. Effects of varying helium abundance 

As mentioned in Sect. ^ M 13 shows strong abun- 
dance variations along its red giant branch, which might 
be attributed to cither deep mixing extending into the 
hydrogen-burning shell (helium mixing) or pollution of the 
stars with the ejecta from an earlier generation of AGB 
stars (helium pollution). Both processes would increase 
the envelope helium abundance of the RGB stars with 
important consequences for the subsequent HB evolution. 
Here we discuss the potential impact of these effects on 
the properties of the blue HB stars in M 13. 



6.1. Helium mixing 

To illustrate the effects of helium mixing, we will use a 
set of mixed sequences computed by Moehler et al. I)2()()()|l 
for a main-sequence helium abundance Y of 0.23 and a 
scaled-solar metallicity [M/H] of —1.54. These sequences 
were evolved up the RGB for different amounts of helium 
mixing using the approach of Sweigart Ijl997a|l . Mass loss 
was included according to the Reimers formulation with 
the mass-loss parameter rjn set equal to 0.45. The evolu- 
tion was then followed through the helium flash to the end 
of the HB phase using standard techniques. The HB locus 
of these helium-mixed sequences in the logg - Toft plane 
is shown in the top panel of Fig. |H1 

Since mixing increases the RGB mass loss due to a 
brighter RGB tip luminosity, a mixed model will arrive 
on the HB at a higher effective temperature than the cor- 
responding canonical model. At the same time mixing in- 
creases the envelope helium abundance in the HB model, 
which leads to a higher energy output of the hydrogen- 
burning shell and hence to a brighter surface luminosity 
(Sweigart & Gross 1976). 

The net effect is to shift the mixed locus in Fig. |H1 to- 
wards lower gravities with increasing effective temperature 
relative to the canonical locus, until a maximum offset is 
reached for 15,000 K < Toff < 20,000 K. At higher tem- 
peratures the mixed locus shifts back towards the canoni- 
cal locus, as the luminosity of the hydrogen-burning shell 
declines due to the decreasing envelope mass. The pre- 
dicted locus along the extreme HB (EHB) does not depend 
strongly on the extent of the mixing, since the luminosities 
and gravities of the EHB stars are primarily determined by 
the mass of the helium core, which is nearly the same for 
the mixed and canonical models. Models at the cool end 
of the mixed locus will also differ little from the canonical 
models, since the cool HB stars will have undergone lit- 
tle mixing on the RGB. This explains why the mixed and 
canonical loci converge at lower temperatures in Fig. |H1 
The size of the offset between the mixed and canonical 
loci depends on the assumed value of rjn, becoming larger 
as r]ji decreases (see Moehler et al. l2000|) . 

6.2. Helium pollution 

The ejecta of AGB stars which have undergone hot bot- 
tom burning should be enriched in helium (Ventura et 
al. l2()()l|l . Thus the pollution of either pre-stellar material 
or existing low-mass stars by this material will increase 
either the overall helium abundance or at least that in 
the envelopes. D'Antona et al. 1)2002(1 have studied the 
evolution of stars that form out of polluted material in 
which the helium abundance was enriched from a canon- 
ical value of r = 0.23 - 0.24 to Y = 0.29. They found 
that while the effects on the pre-HB evolution are quite 
small, the HB morphology itself is strongly affected by he- 
lium pollution. The main effect is to shift the HB stars to 
hotter temperatures due to smaller turnoff masses for a 
given cluster age and to increased mass loss on the RGB. 
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Fig. 8. Temperatures and gravities of the programme 
stars in M 3 and M 13 (from metal-rich model atmo- 
spheres for stars hotter than 12,000 K and metal-poor 
model atmospheres for cooler stars, in both cases includ- 
ing metal lines) compared to non-canonical evolutionary 
models. The zero-age and terminal-age HB predicted by 
canonical models for Y = 0.23 and [M/H] = -1.54 are plot- 
ted as solid lines. 

a) The zero-age and terminal- age HB for [M/H] = —1.54 
computed with helium mixing for a Reimers mass-loss pa- 
rameter riR = 0.45 (see Sect. IG.ll for details). 

b) The zero-age and terminal-age HB computed with he- 
lium pollution for helium abundances of y = 0.28 and 
Y = 0.33 (see Sect. lO for details). 



Polluted stars should thus be found primarily among the 
hottest HB stars. Indeed, D'Antona et al. (2002,i) suggest 
that helium pollution is responsible for the extended blue 
tails found in such clusters as M 13 and NGC 6752. 

To investigate the effects of helium pollution, we com- 
puted two additional sets of HB sequences for [M/H] ~ 
— 1.54 in which the helium abundance was increased by 
0.05 and 0.10 above our canonical value of y = 0.23. The 
HB loci defined by these helium-polluted sequences are 
plotted in the bottom panel of Fig. |S1 The gravity offset 
between these loci and the canonical locus is primarily due 
to the luminosity difference between the helium-polluted 
and canonical models. 



HB stars derive their energy from two sources: he- 
lium burning in the core and hydrogen burning in a shell. 
Since an increase in the initial helium abundance leads 
to a decrease in the helium-core mass (Sweigart & Gross 
11978(1 . the energy output of the helium core will be lower 
in a helium-polluted model than in a canonical model. 
Along the EHB, where hydrogen-shell burning is unim- 
portant, a helium-polluted star will therefore be fainter 
and have a higher gravity than a canonical star. The op- 
posite is true at cooler temperatures where the hydrogen- 
burning shell is a major energy source. The higher en- 
velope helium abundance of a helium-polluted star then 
leads to a brighter surface luminosity and thus to a lower 
gravity. At the transition between the blue HB and the 
EHB (Toff « 20,000 K), the gravities of the canonical and 
helium-polluted models are virtually identical. This is the 
point where the higher hydrogen-burning luminosity of a 
helium-polluted star offsets the lower helium-burning lu- 
minosity. 

One aspect of the helium-polluted loci in Fig. [S] 
requires further comment. According to the helium- 
pollution scenario outlined by D'Antona et al. H2002|l . the 
stars at the cool end of an M 13-like HB, i.e., those near 
the top of the blue tail, would not be helium polluted. 
Thus in an actual cluster one would expect the cooler 
HB stars to lie within the canonical locus. However, as 
one goes to higher temperatures along the blue tail, the 
fraction of the HB stars that are helium polluted would 
increase, and thus the locus predicted by the D'Antona 
et al. scenario would shift away from the canonical locus 
towards the lower gravities of the helium-polluted loci in 
Fig. IHl This shift to lower gravities is not seen in Fig. |H1 
because the helium-polluted loci in this figure assume that 
all of the HB stars including the cooler stars are helium 
polluted. 



6.3. Comparison with observations 

In Figs. IHl and El we compare the effective temperatures, 
surface gravities and masses derived from metal-rich (Tea 
> 12,000 K) and metal-poor (Tch < 12,000 K) model at- 
mospheres to the scenarios described above. While the 
helium-mixed models provide a better description of the 
temperatures and gravities derived for the stars between 
12,000 K and 16,000 K than the canonical models, they do 
not resolve the problem with the low observed masses. The 
same is also true for the helium-polluted models, where we 
again find that the observed masses are systematically too 
low. The fact that neither the helium-mixed nor helium- 
polluted models can account for the low masses raises 
doubt about the reality of any gravity and/or tempera- 
ture discrepancies between the programme stars in Fig.|Sl 
and canonical predictions. Without a better understand- 
ing of possible systematic errors, it is not possible to argue 
for or against the above noncanonical scenarios, given the 
present data. 
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Fig. 9. Temperatures and masses of the programme stars 
in M 3 and M 13 (from metal-rich model atmospheres for 
stars hotter than 12,000 K and metal-poor model atmo- 
spheres for cooler stars; in both cases metal lines were 
included in the theoretical spectra) compared to evolu- 
tionary tracks. 

a) The solid line marks the canonical ZAHB and the 
mixed ZAHB is given by the short-dashed line (see 
Sect. IB.ll for details). 

b) Again the solid line marks the canonical ZAHB and 
the polluted ZAHBs are given by the dotted (Y=0.28) 
and long-dashed Hne (Y=0.33), respectively (see Sect. 16.21 
for details). 

7. Conclusions 

We have obtained low-resolution (3.4A) spectroscopy of 
22 hot HB candidates in M13, and four in M3, in order 
to derive atmospheric parameters (effective temperatures 
and surface gravities) as well as abundances of helium, 
magnesium, and iron. One star (G43) in M13 turned out 
to be a UV-bright star, while the remaining targets appear 
to be bona-fide HB stars. For the stars between 10,000 K 
and 12,000 K, the atmospheric parameters derived from 
fitting the Balmcr lines and from Stromgren photometry 
are in good agreement. For the stars cooler than 10,000 K, 
there is a discrepancy between parameters derived from 
Stromgren photometry and those derived from fitting the 
Balmer lines. We suspect that there are problems in the 
Balmer line fitting because the derived temperatures show 



a gap between about 9000 K to 10,000 K, and because the 
derived iron abundances show a wide scatter. However, 
despite numerous tests, we were unable to determine the 
cause of the disagreement in the derived parameters. 

For stars hotter than 12,000 K in both clusters, we 
find evidence for helium depletion and a large iron enrich- 
ment, consistent with the results on M13 by Behr et al. 
ifTM^ and on NGC 6752 by Moehler et al. PTOjl . The 
similar temperatures for the onset of radiative levitation in 
M3, M13, and NGC 6752 suggest that this phenomenon is 
unrelated to the HB morphology. Instead, the onset of ra- 
diative levitation may be connected to the disappearance 
of the surface convection zone, as suggested by Sweigart 
(2002). 

We compare our temperature and gravity results with 
the predictions of the helium mixing scenario of Sweigart 
l|1997b|l and the helium pollution scenario of D'Antona 
et al. H2002|l . These scenarios are attractive because they 
can explain the origins of the HB blue tail in globular 
clusters such as M13 without requiring a fine-tuning of 
the mass loss, and because they can relate the blue tail 
to the observed abundance anomalies on the RGB. Both 
scenarios predict lower gravities (and larger luminosities) 
for stars near 15,000 K. For stars cooler than 12,000 K, 
we find the observed gravities in agreement with canonical 
models. This result is consistent with the work of Caloi 
1)200 ip . who compared the HB luminosities of M3 and M13 
to conclude that there was no evidence for a substantial 
surface helium abundance increase for HB stars near the 
temperature of the RR Lyrae stars. 

For stars hotter than the onset of radiative levita- 
tion at 12,000 K, we fit the Balmer lines using metal-rich 
([Fe/H] — -1-0.5) model atmospheres to derive the temper- 
ature and gravity. Although the use of metal-rich atmo- 
spheres reduces the discrepancy with canonical models, we 
still find an offset to lower gravity of about 0.2 dex for stars 
between 12,000 K and 16,000 K, similar to what was seen 
at temperatures hotter than 15,000 K in NGC 6752 by 
Moehler et al. f 2000') . However, there are a couple of rea- 
sons why this result should be viewed with caution. First, 
the derived masses for these stars appear to be too low 
compared to either the canonical or the helium-enriched 
scenarios. Second, the metal-enriched atmospheres used 
to derive these parameters are not well-tested, and do 
not take into account the strongly non-solar abundance 
ratios. Future high-resolution spectroscopy could help re- 
duce these possible systematic errors in gravity, by allow- 
ing the atmospheric parameters and abundances to be de- 
termined iteratively, using model atmospheres with non- 
solar abundance ratios. 
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Fig. 2. Normalized spectra of the programme stars in M 13 and M 3. The part shortward of 3900 A was normalized 
by taking the highest fliix point as continuum value. 



